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Abstract. While real Hamiltonian mechanics and Hermitian quantum mechanics 
can both be cast in the framework of complex canonical equations, their complex 
generalisations have hitherto remained tangential. In this paper quaternionic and 
coquaternionic (split-signature analogue of quaternions) extensions of Hamiltonian 
mechanics are introduced, and are shown to offer a unifying framework for complexified 
classical and quantum mechanics. In particular, quantum theories characterised by 
complex Hamiltonians invariant under space-time reflection are shown to be equivalent 
to certain coquaternionic extensions of Hermitian quantum theories. One of the 
interesting consequences is that the space-time dimension of these systems is six, not 
four, on account of the structures of coquaternionic quantum mechanics. 

PACS numbers: 03.65.Aa, 02.30.Fn, 03.65. Ca 

This paper concerns the relation between complexified classical and non- Hermitian 
quantum mechanics, and their surprising hnks to quaternionic and coquaternionic 
mechanics. The main finding is that complexified mechanical systems with real energies 
studied extensively in the hterature over the past decade can alternatively be thought 
of as certain coquaternionic extensions of the underlying real mechanical systems. This 
identification leads to the possibility of employing algebraic techniques of quaternions 
and coquaternions to tackle some of the challenging open issues in complexified classical 
and quantum mechanics. 

Complex (i.e. non-Hermitian) Hamiltonians have long been employed to describe 
open quantum systems, decay and scattering phenomena [Tj. Further, since the 
realisation that complex operators respecting space-time refiection (PT) symmetry may 
possess entirely real spectra [2], there have been considerable research interests in 
examining both physical and mathematical properties of quantum systems described 
by non-Hermitian Hamiltonians with real spectra. More recently, the interest in 
these systems has increased notably, in part owing to experimental realisations of the 
phenomenon of the PT phase transition and other theoretically predicted effects |3l H] . 

Complexified classical mechanics has also been studied intensely both in the context 
of semiclassical calculations and as a classical analogue of non-Hermitian quantum 
mechanics [5], El [7]. At classical level, complex-extended mechanics is characterised 
by a phase space spanned by complex canonical variables. The quantum analogue of a 
phase space — in the strict sense of a symplectic manifold upon which the dynamics is 
governed by Hamilton's equations of motion — is the space of pure states, i.e. the space 
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of rays through the origin of the Hilbert space. On this ray space, the wave equation of 
Schrodinger in Hermitian quantum mechanics can be expressed as Hamihon's equations 
of motion [HI |9]. An apparent puzzle in this context is that the analysis of complex 
extensions of quantum mechanics has thus far been confined to the real state space 
fiO[ im [12] (with the exception of ^3]). Therefore, the nature of complexification 
considered so far for quantum mechanics has a rather different characteristic from that 
for classical mechanics. If classical phase spaces can be extended into the complex 
domain, it seems paradoxical that quantum state space cannot be extended in an 
analogous manner. (Quaternionic quantum mechanics [TU [15] can be viewed as 
representing 'complex' extensions of the underlying real mechanics for real (Hermitian) 
quantum systems, although this is not the point of view commonly adopted.) 

In the present paper we propose a new approach to these apparently different 
complex generalisations by making use of quaternionic and coquaternionic formulations. 
The key idea is that phase-space variables consist of canonical conjugate pairs; hence a 
complexified mechanical system necessarily involves pairs of complex numbers. On the 
other hand, pairs of complex numbers can usefully be treated as a single quaternion. 
It is natural therefore to enquire whether complexified mechanics can be represented 
more concisely in terms of quaternions, inasmuch as real phase space variables can 
more conveniently be represented as a single complex number. A deeper motivation for 
such a study arises from the fact that symmetry properties of quaternions are closely 
related to symmetries of (Euclidean) space-time [16]; while complexified dynamical 
systems discussed here are characterised by Hamiltonians invariant under space-time 
reflection (PT) symmetry. The appearance of space-time symmetry in both complexified 
mechanics and quaternionic algebra suggests that the latter might constitute a useful 
tool for the investigation of the former. A further motivation comes from the observation 
that if the parity operator is assumed trace free, then the real dimensionality of the 
Hilbert space of a PT-symmetric quantum system is a multiple of four, not two |17j . 
This suggests that quaternions a priori are suitable candidates for the characterisation 
of phase space variables for such theories. 

With these observations in mind, we shall explore the possibility that complexified 
mechanics can alternatively be viewed as a version of quaternionic mechanics. We 
shall find, perhaps somewhat unexpectedly, that PT-symmetric systems are in fact 
related to coquaternions, rather than quaternions. In this connection it is worth 
remarking that symmetries of coquaternions are related to the Lorentz group, rather 
than Euclidean group, in the sense that every rotation in the Minkowski three-space 
can be expressed in terms of coquaternions. Hence PT symmetry a priori has more 
in common with Lorentzian space-time than Euclidean space-time. It follows that PT- 
symmetric quantum mechanics, although related, are not equivalent to the traditional 
quaternionic quantum theories of [TU [15]. 

Our hope is that the approach introduced here may ultimately help to tackle some 
open issues in the area of complexified mechanics that have only been answered at 
most partially. These include, among others: The classical-quantum correspondence 
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of complexified mechanics [TSl dSl 1201 EI]; tfie derivation of statistical mechanics of 
classical PT-symmetric systems [22]; the characterisation of combined systems (e.g., 
entanglement) in PT-symmetric quantum mechanics, and the formulation of a theory 
of space-time that is consistent with complex Hamiltonians. 

Quaternions and coquaternions. Although properties of quaternions are well known, 
we find it convenient to briefly summarise here some of the key relations. The algebra 
of quaternions involves three imaginary units, denoted here by j, /c, satisfying the 
relation 

i'=3^ = k^ = i3k = -l, (1) 

carved into the stone of the Brougham Bridge by Hamilton in an act of "mathematical 
vandalism" [23j. It follows that quaternions satisfy the following cyclic relation: 

ij = —ji = k, jk = —kj = i, ki = —ik = j. (2) 

Hence matrix representations of the three imaginary units are just the Pauli matrices 
multiplied by the complex number. 

A generic quaternion can be represented in the form q = qo + iqi + jq2 + kqs, where 
{qi}i=o,i,2,3 are all real. It is not difficult to see that quaternions that commute with all 
other quaternions are reals, and that the totality of quaternions that commute with a 
given quaternion forms a subset isomorphic to complex numbers. The conjugate of a 
quaternion q is given by g = go ~ ^Qi ~ 312 — kq^- It follows that the squared modulus 
of g is gg = go + 9'i + ?2 + ^'35 that the inverse of g is g^^ = g/(gg). A quaternion 
can be expressed in polar form 

q= \q\e^i^i = |g| (cos 6*5 + sin 615), (3) 

where iq = {iqi + jq2 + kq3)/\q\ and 9g = cos"^(go/|g|). 

The symmetry properties of quaternions are thus closely related to the group SU (2). 
The SU{1, 1) analogues of quaternions are the coquaternions [2l|, also known as split 
quaternions. They satisfy the relation 

= -1, f = e = ijk = +1 (4) 

in place of ([T]), and the skew-cyclic relation 

ij = —ji = k, jk = —kj = —i, ki = —ik = j (5) 

in place of ([2]). Like quaternions, the conjugate of a coquaternion g = go+^gi+jg2+^Q'3 is 
g = go— ^gi— jg2— fcga- It follows that the squared modulus of a coquaternion is indefinite: 
= Qo ~^ ~ (I2 ~ Qs- Unlike quaternions, a coquaternion need not have an inverse 
g~^ = q/{qq) if it is null, i.e. if gg = 0. The polar decomposition of a coquaternion thus 
involves trigonometric and hyperbolic functions in the usual way. Unlike quaternions, 
however, coquaternions do not admit finite-dimensional unitary representations. 

Complex formulation of real classical and quantum dynamics. Before we proceed 
to explore the quaternionic or coquaternionic extensions, it is useful to represent the 
dynamical theory of elementary mechanical systems in terms of complex phase-space 
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variables We start with a one-dimensional dynamical system, described 

by real canonical variables p and x, and introduce the complex conjugate variables 
z = (x + ip) / \/2 and z = {x — ip)/ \/2. The canonical equation of motion then reads 
dz OH 

'Tt = 9F' 

where H = H{z, z) denotes the Hamiltonian. In this formulation it becomes apparent 
in which way quantum Schrodinger dynamics can be viewed as a special case of 
Hamiltonian mechanics [25]. Consider a quantum system characterised more generally 
by a state vector 1^;), with components z„. The Schrodinger equation, with a Hermitian 
Hamiltonian H, can then be written in the form izn = dH/dzn, where H = {z\H\z) = 
J2mn HmnZnZm IS the expectation value of H. In particular, the Hamiltonian function 
of a quantum-mechanical system is bilinear in the canonical variables, that is, it is 
essentially the Hamiltonian of coupled harmonic oscillators. 

With this in mind, we shall devote a great part of the present paper to the 
investigation of the bilinear Hamiltonian H = zz, leading to the equation of motion 
idz/dt = z. The trajectories in the z-plane are concentric circles about the origin, 
traversed with a constant angular velocity, which define a linear critical point called 
a centre. If we perform a Wick rotation, we obtain dz/dt = z (or dz/dt = — z), 
which defines a linear critical point called a focus. The trajectories are rays leaving 
(or approaching) the origin, traversed with radially increasing velocity. More generally, 
consider the equation dz/dt = bz, where 6 is a complex constant. If both the real and 
imaginary parts of b are nonzero, then the corresponding linear critical point is a vortex, 
and the orbits consist of concentric spirals leaving (or approaching) the origin traversed 
with equal radially increasing velocities. Quantum mechanically, all fixed points of 
unitary dynamics are centers, whereas in a PT-symmetric quantum theory a pair of 
centres can turn into a pair of foci or vortices upon symmetry breaking. 

Complexified classical and quantum mechanics. Let us now summarise briefly 
the conventional way in which classical and quantum mechanics are extended into 
the complex domain. The complexified classical mechanics typically starts from the 
real canonical equations of motion p = —dH/dx and x = OH /dp, and considers 
the complexification p = po + ipi and x = xq + ixi. In terms of the four real 
variables {pQ,Xo,pi,Xi) the Hamiltonian now takes the form H = Ho{pq,xq,Pi,Xi) + 
iHi{po,xo,pi,Xi). Assuming that H is analytic and thus satisfies the Cauchy-Riemann 
conditions dHo/dxo = dHi/dxi, dHo/dxi = —dHi/dxo, dHo/dpo = dHi/dpi, and 
dHo/dpi = —dHi/dpo, the equations of motion reduce to 

dHo . dHo . dHo . . dHo 

Po = — , xo = - — , Pi = - — , and xi = (7) 
OXq opo OXi opi 

In particular, since the total energy Hq + iHi is conserved under the Hamiltonian 

dynamics, if the initial energy is real so that Hi = 0, then it will remain real. Thus, 

a complexified Hamiltonian system is equivalent to a real two-dimensional system for 

which the position and momentum are interchanged in the second dimension [5l E] . 
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The investigations of 'complexified' quantum mechanics, on the other hand, 
consider complex Hamiltonians of the form H = Hq + iHi, while typically leaving 
the quantum phase space real. The equation of motion is still given by the Schrodinger 
equation. That is, writing Zn = {xn + ipn)/V^ for the components of the state vector, 
we have izn = dH/dzn, where the Hamiltonian function is the expectation value of H, 
which is now complex. It is then straightforward to verify that equations of motion for 
the real phase space variables Pn, Xn are given by 

dHo dHi , . dHo dHi 
Pn = ^ — and Xn = ^ ^ — . (8) 

OXn OPn OPn OXn 

This structure is a combination of a Hamiltonian symplectic flow generated by Hq and a 
Hamiltonian gradient flow generated by Hi, and also appears in the semiclassical limit 
of non-Hermitian quantum dynamics |18j . 

The dynamics of complexified quantum systems considered in the literature, for 
the most part, are characterised by ([H]), and not by (^^. This is the sense in which 
the nature of complexification in classical mechanics has been different from that of 
quantum mechanics. To see in which way these two formalisms might be unified by a 
quaternionic or coquaternionic approach, we shall now investigate the special case of a 
bilinear Hamiltonian, before turning to the more general case. 

Quaternionic and coquaternionic oscillators. We again examine a one- dimensional 
dynamical system with the Hamiltonian H = zz, but now z is assumed a quaternion. 
That is, although the dynamical variable is given hy z = {x + ip) / a/2, we let p and 
X extend into the complex domain according to the prescription p — > Po + jPi aiid 
X — )■ XQ+jxi. It follows from the algebraic property ([2]) that z = (xo+ipo+jxi+kpi) / \/2. 
We assume for simplicity a superselection rule [2] that singles out the complex number 
i for determining the direction of time. Hence the canonical equation of motion is given 
by (E]) with H = pI + pf + xl + xf, and we obtain 

Po = -xo pi = -xi 

Xo=Po Xi=pi. 

Thus, a quaternionic oscillator is equivalent to a two-dimensional oscillator. We remark 
that working with the real variables {p,x), a natural way of extending the theory into 
quaternionic domain would not have been apparent. 

Let us compare this system with a complex PT-symmetric oscillator [7], where we 
take the Hamiltonian H = + x"^) for which equations of motion read p = —x and 
X = p, and extend the phase space variables into the complex domain: p — )■ po + Wi and 
X — )■ xo + ixi. The resulting dynamical equations are given by ([9]). While equations of 
motion agree, and indeed they share the same fixed-point structure, these two theories 
are nevertheless distinct because of boundary conditions: The energy E = Pq+pI+Xq+xI 
is positive definite for quaternionic oscillators; whereas the energy E = pl—pl + x^—xl 
is indefinite for PT-symmetric oscillators. 

Next we consider a coquaternionic oscillator, with the same superselection rule. (It 
is worth remarking that unlike quaternionic quantum mechanics where one has to make 
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the choice for the direction of time, in coquaternionic quantum mechanics the imaginary 
unit i is naturally preferred on account of the fact that P = —1, while = fc^ = +1.) 
In this case, equations of motion remain the same, but owing to the split-signature 
of coquaternions the energy is given hj E = pi — pf + — xl- Hence we conclude 
that a complex PT-symmetric oscillator is equivalent to a coquaternionic oscillator. An 
important point to note here is the fact that in the complex formulation of mechanics 
the construction of the Hamiltonian such as H = zz implicitly involves the notion of an 
inner product, which is not evident in Hamilton's formulation involving real variables 
{p,x). In the case of a complex PT-symmetric oscillator, this inner product is in effect 
determined by a PT conjugation, which, owing to the fact that the parity operator is 
trace free, has a split signature equivalent to a coquaternionic norm. 

At the level of classical mechanics, the appearance of a PT inner product seems 
unproblematic because it has no impact on equations of motion. However, if we attempt 
to formulate statistical mechanics, for instance, then the use of a PT inner product 
raises a severe obstacle: For any bounded energy, the energy shell in phase space is not 
compact, hence an equilibrium microcanonical distribution cannot be defined. At the 
quantum level, the issue of indefiniteness in a PT inner product has been recognised 
earlier because inner products are related to probabilities, which have to be nonnegative. 
To remedy this issue, an alternative inner product, based on a CPT conjugation, has 
been introduced [271 EH]- One possible way forward in the case of a complex classical 
oscillator therefore is to modify the inner product such that indefinite components pick 
up an additional minus one. Then energy shells for finite energies become compact, 
allowing, in particular, for a rigorous formulation of statistical mechanics. 

Coquaternionic mechanics with real and complex energies. More generally, consider 
now a real-valued Hamiltonian H{z,z) where z = {xq + ipo + jxi + kpi)/\/2 is a 
coquaternion, and hence \/2d/dz = d/dxo — id/dpo + jd/dxi + kd/dpi. From ([6]) 
we thus deduce that 



zz 

and hence that 



dH .dH .OH ^dH 

dxo dpo ^ dxi dpi 



(10) 



dH dH dH ^ dH 

Po = = — , pi = — , and Xi = (11) 

dxQ dpo dxi dpi 

Comparing these equations with ([7]) where the initial energy is assumed real, we 
find therefore that a PT-symmetric classical system with real energy is equivalent 
to a coquaternionic mechanical system with real energy. Note that in a quaternionic 
mechanical system the role of pi and Xi are interchanged, and thus it is merely equivalent 
to a two-dimensional real system. 

An alternative case to consider is where the phase space or the state space variables 
are kept real so that z remains a complex number, but the energy H = Hq + iHi is 
made complex. In this case, it follows from (|6]) that equations of motion read 

dH, dH^ . dH.dH^ 

P = — 7, h -7T— and X = — h (12) 

dx dp dp dx 
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which is of course a special case of ([S]). If we now allow both energy and phase space 
variables be complex, then in the coquaternionic case, we obtain: 

dHo dHi dHo dHi 
Po = ^ ' ^0 = ^ h ^ — , 

OXq Opo Opo OXq 

OXi Opi Opi OXi 

These equations constitute natural generalisations of f|T2l) in the sense that they embody 
the structure of a combination of a Hamiltonian symplectic flow generated by Ho and a 
Hamiltonian gradient flow generated by Hi. 

We remark, more generally, that if the Hamiltonian is coquaternionic-valued so that 
H = Ho + iHi+ jH2 + kH3, then ([10]) implies 

dHo dHi dH2 dH^ . dHo dHi dH2 dH^ 
Po = — t: \- ^ \- ^ \- ^ — , xo = ^ h 



dxo dpo dxi dpi dpo dxo dpi dxi 



. dHo , dHi dH2 , dH^ . dHo , dHi dH^ , dH^ 
^^=9^7+9^-9^^+9^' "^ = -9^+9^7+9^+9^^ (''^ 
for the dynamical equations. 

Two oscillators and two-level systems. The orbit space of a classical system of a pair 
of oscillators is equivalent mathematically (though not physically) to a quantum two- 
level system. In particular, the reduced phase space of both systems is just the Bloch 
sphere S"^. This can be seen by noting that the Hamiltonian of a pair of uncoupled 
oscillators in complex coordinates is H = z^z^ + ^2^2- If the energy is fixed we obtain 
the three-sphere ^12:1 + Z2Z2 = 1, upon which the Hamiltonian flow acts as scalar 
multiplication Zn — e**2;„ [29|. Hence the orbit space is just the two-sphere resulting 
from the Hopf map S*^ — )■ S*^. 

Similarly, a pair of classical quaternionic oscillators is equivalent to a quaternionic 
two-level quantum system. In this case, once energy is fixed we obtain the seven- 
sphere ziZi + ^2-2:2 = 1, upon which the Hamiltonian flow acts as quaternionic scalar 
multiplication. Hence the associated phase space (Bloch sphere) is a four-sphere S'^ 
resulting from the Hopf map 5*^ — )■ 5"^. To see the structure of the state space in the 
quaternionic case more explicitly, it suffices to note that a generic normalised state 
vector in a quaternionic two-level system can be written in the parametric form 

where = i cos 02 +j sin 02 cos (p^ + k sin 02 sin 03 . It is then straightforward to see that 
(^, 01, 02, 03) constitute spherical coordinates for an 5^ in R^. In higher dimensions, the 
state space is just the quaternionic projective Hilbert space endowed with a quaternionic 
Fubini-Study metric arising from transition probabilities. In the case of a pair of 
coquaternionic oscillators, the resulting orbit space is hyperbolic (a four-dimensional 
analogue of the Poincare disk, which is an example of a Siegel domain |30]). 

In the case of a quaternionic or coquaternionic two-level system, a generic Hermitian 
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Hamiltonian can be expressed in the form 

5 

H = liol + J2^i^h (16) 

1=1 

where a; e R, n is a unit vector on 5'^ C R^, and 




0"! 



(17) 

are the (co)quaternionic Pauh matrices. For a quaternionic Hermitian matrix the 
eigenvalues are real; whereas they are either real or appear as complex conjugate pairs 
in the case of a coquaternionic Hermitian matrix. In either case, for two-level systems 
there are six parametric degrees of freedom (consistent with the observation that the 
most general PT-symmetric 2x2 Hamiltonian has six exogenous parameters [31]). In 
standard quantum mechanics, it follows from the construction of spin-orbit interaction 
(the so-called Pauli correspondence) that one can identify the spin operators with the 
three spatial directions. Similarly, (fT7|) suggests that in a quaternionic or coquaternionic 
quantum mechanics, the ambient space has five dimensions. 

It is worth remarking that owing to the lack of commutativity of the imaginary units 
there are "observable" effects in quaternionic or coquaternionic mechanics that have no 
analogue in real mechanics. For instance, a Hamiltonian proportional to an identity 
matrix can generate nontrivial dynamics [32]. In particular, if we let i determine the 
preferred complex structure, and take the Hamiltonian H = then the dynamics 

gives rise to a Rabi oscillation between the (J4 and (T5 directions. This can be seen most 
easily in the Heisenberg picture: 

cf^it) = cos((X't)(T4 + sin (wt) 0-5, (18) 

where we have chosen the initial state to be the (T4 eigenstate. Therefore, we can generate 
dynamics with vanishing energy gap. This is reminiscent of the so-called "arbitrary fast" 
quantum state transport effect [TUl [HI [12] . 

Let us now consider a 2 x 2 coquaternionic Hermitian Hamiltonian and express this 
in the form 

^=1^"^ VI, (19) 




where q = qo + iqi + jq2 + kq^. The eigenvalues are E± = s ± \/t^~+~^, where qq = 
go + 9i — Q'i ~ Q's • Therefore, in the coquaternionic representation, if t'^ + qo + qf > + li^ 
then the eigenvalues are real, which can be identified as the region of unbroken PT 
symmetry; otherwise, the eigenvalues form a complex conjugate pair. 

Working in the basis for which H is diagonal, the phase-space function for the 
Hamiltonian is if = Z1E+Z1 + Z2E-Z2 (note that in general E± and Zi need not commute. 
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since Ej- need not be real). It follows that equations of motion are given by ii = iE^zi 
and Z2 = iE_Z2- In particular, if PT symmetry is unbroken, then eigenvalues are real 
and the fixed points associated with the dynamical evolution are centres; whereas if the 
symmetry is broken, then the fixed points form a vortex pair. This is a characteristic 
feature of the so-called PT phase transition whereby a pair of centres coalesce and turn 
into a pair of vortices (a sink and a source) |33]. Such a transition does not occur in 
the case of a quaternionic quantum theory, and this establishes the inequivalence of 
PT-symmetric and quaternionic quantum theories. 

Discussion. We have shown that PT-symmetric quantum theories are equivalent 
to certain coquaternionic extensions of Hermitian quantum theories, at least in finite- 
dimensional cases. This equivalence suggests that the underlying space-time dimension 
is six, not four, on account of the two additional dimensions generated by the additional 
imaginary units. This follows from the fact that spin-orbit interaction can only be 
defined for five spatial dimensions in quaternionic or coquaternionic quantum theories. 

Needless to say, there is a substantial hterature dealing with various aspects 
of six-dimensional space-times. Apart from the practicality that calculus on a six- 
dimensional space facilitates various analysis in four dimensions |3ll [351 ES] , models for 
six- dimensional space-times have also been proposed at a fundamental level. Empirical 
evidences for extra dimensions are often sought at inaccessibly high energies. However, 
by identifying physical phenomena predicted by coquaternionic quantum theories that 
cannot be described by standard complex quantum mechanics (such as the Rabi 
oscillation with vanishing energy gap described above), it is possible to seek compelling 
evidences for the existence of extra dimensions at low energies (cf. [37] for a discussion 
on difficulties in detecting quaternionic manifestations in scattering experiments). To 
this end, it is worth remarking that evidences for octonionic quantum mechanics would 
suggest that the underlying space-time has dimension ten. 

Some of the challenging open issues of PT-symmetric quantum theories include 
the formulation of combined systems to characterise such notions as entanglement 
and coherent states. The mathematical difficulty might be rooted in the fact that 
quaternionic projective Hilbert spaces (unlike complex projective space of ordinary 
quantum mechanics [38]) do not admit algebraic varieties that characterise subspaces 
of particular physical characteristics [39]. Nevertheless, progress has been made in 
the context of quaternionic quantum mechanics to formulate combined systems [IQ] or 
coherent states [HI Hlj . It may be that these techniques can be applied to investigate 
open issues associated with combined systems for PT-symmetric quantum theories. 

EMG is supported by an Imperial College Junior Research Fellowship. 
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